Field induced tricritical behavior in the S* = 1/2 quasi one-dimensional frustrated Ising 
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The results of extensive histogram cluster heat-bath Monte Carlo simulations on the critical 
behavior of the quasi-one dimensional Ising antiferromagnet on a stacked triangular lattice are 
presented. A small applied field is shown to induce a crossover from XY universality to mean- 
field tricritical behavior. Experimental estimates of critical exponents suggest that these two types 
of phase transitions are observed in S = 1 CsNiCla and 5=1/2 CsCoBra, respectively. The 
present results demonstrate that this difference can be explained by an unusual staggered magnetic 
field arising from quantum exchange mixing previously proposed to account for spin excitations in 
5=1/2 quasi-one-dimensional Ising antiferromagnets. 
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Meaningful comparisons between experimental and 
theoretical results for critical exponents can be notori- 
ously difficult due to inadequacy of standard techniques 
to sample thermodynamic quantities close enough to the 
phase transition region. It is often the case that small 
perturbations such as weak anisotropy govern the uni- 
versality of critical phenomena but only very close to 
the critical temperature (T^r), rendering access through 
simulations or measurements problematic. Such issues 
are exacerbated by critical behavior thought to be influ- 
enced by magnetic frustration pL]. Additional complica- 
tions can arise in model systems where exchange interac- 
tions are strongly anisotropic or in cases where quantum 
effects govern spin dynamics. The Ising antiferromagnet 
on a stacked triangular lattice (ISTAF) with strong c- 
axis exchange and weak inter-chain exchange interactions 
is an example of such a model system. Many quasi-one- 
dimensional (ID) ABX3 compounds exhibit strong short- 
range intra-chain magnetic order (SRO) at temperatures 
well above the onset of true three-dimensional (3D) long 
range order (LRO) 0. The present work demonstrates 
that quantum effects which govern ID spin dynamics can 
impact effective critical properties of the 3D system. 

Exchange mixing due to an unquenched orbital mo- 
ment and crystal field effects in quasi-lD S = 1/2 
ISTAF's has been proposed to give rise to an unusual 
contribution to the effective Hamiltonian in the form of 
a staggered magnetic field with a periodicity of two along 
the c-axis. This term has been employed to explain low 
temperature neutron scattering spectra in CSC0CI3 and 
CsCoBrg fs^l and Raman data in TIC0CI3 1- Invoking 
such an effect relies on the well-developed SRO intrinsic 
to quasi-lD magnetic systems at low T. Its applicabil- 
ity has been questioned, however, due to the fact that 
the inclusion of such a term in the model Hamiltonian 
essentially assumes ID LRO and is thus a mean-field ap- 
proximation inappropriate for soliton dynamics fH]. The 
impact of this field term on the measurements of effective 



critical behavior observed close to Tn is not subject to 
this restriction and deserves consideration even though 
it may be strictly irrelevant in the renormalization group 
sense at T^. In this work, the effect of this quantum 
staggered field in the S — 1/2 quasi-lD ISTAF is mim- 
icked in Monte Carlo (MC) simulations by considering 
the influence of a small uniform magnetic field applied 
to a model with ferromagnetic intra-chain interactions 
with AF inter-chain coupling. No such effective quan- 
tum staggered field is predicted to occur in the case of 
non-5 = 1/2 STAF's such as (5 = 1) CsNiCls which has 
been shown to exhibit integer-spin Haldane-gap phenom- 
ena Q. 

Critical behavior of the ISTAF have been investigated 
using numerous theoretical approaches, MC simulations 
and experimentally. Despite the large number of stud- 
ies there is still no clear consensus in the literature re- 
garding the critical exponents characterizing the phase 
transition between the partially disordered phase (PD) 
where one of the three sublattices remains disordered, 
and the paramagnetic phase, at T^. Most modeling re- 
sults [3, ^ support the notion from symmetry arguments 
that the transition at T/v belongs to the XY universal- 
ity class [8]. There have been suggestions, however, from 
several MC simulations that critical exponents are close 
to mean-field tricritcal [l^] , but limited statistics and 
data analysis were used in these investigations. From the 
experimental standpoint, a number of past neutron scat- 
tering studies of CsCoBrs and CsCoCls [2] have obtained 
sets of critical exponents that largely support the idea of 
XY critical behavior. Intriguingly, recent high resolu- 
tion neutron scattering experiments on CsCoBrs by Mao 
et al. IlhI revealed results which suggest with tricritical 
mean-field behavior and it was speculated that this could 
be attributed to the anisotropic nature of the exchange 
interaction or a consequence of the quantum nature of 
the S — 1/2 spins. Experimental data on weakly axial 
S = 1 CsNiCla with a quenched orbital moment Q are 
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consistent with XY universality 0, Both types of 
experimental systems, strongly Ising and weakly Ising, 
exhibit magnetic transitions at Tjv with the same sym- 
metry and thus should belong to the same universality 
class. 

Motivated by the new experimental data of Mao et 
ai, extensive MC simulations using the efficient Cluster 
Heat Bath algorithm (CHB)[i3, |3, combined with the 
histogram methodflB]. were used here to investigate the 
effects of anisotropic exchange as well as a small applied 
field on the critical properties of the classical ISTAF. 
This work serves to the extend previous MC simulations 
of the ISTAF mentioned above, but especially that of 
Koseki et al. |10|] , where the CHB was used with strongly 
anisotropic exchange but with no applied field, as well as 
that of Netz and Berker [l^ and Plumer and Mailhot 
01 who considered isotropic exchange interactions with 
a nonzero applied field. The present results also compli- 
ment MC simulations on the frustrated XY-STAF where 
anisotropy in the exchange interactions was shown to in- 
duce a first-order transition in agreement with recent 
experimental data 

MC simulations in the present work used the order- 
parameter culumant crossing method[l3] to determine 
the critical temperature of the system. Finite-size scal- 
ing analysis at T/v was then used to extract estimates 
of the critical exponents and i'. Both of these tech- 
niques have been thoroughly tested on frustrated Ising 
and Heisenberg spin models [l^ . 

For the purposes of determining the classical critical 
behavior it is adequate to use a model system with ferro- 
magnetic exchange interactions along the c-axis, Jq , with 
or without a small applied uniform field H directed along 
the c-axis. In addition to near-neighbor AF interactions 
between chains, Ji, the effects of a small next- nearest- 
neighbor exchange coupling, J2, was also considered as 
in previous MC simulations [21|. These effects can be 
incorporated in the Ising Hamiltonian: 



(1) 
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where the spin at any site is = ±1. In zero applied field 
4 different cases were investigated with parameters listed 
in table HI For all cases Jq and J2 are ferromagnetic and 
Ji is antiferromagnetic. Finite-size scaling analysis was 
also performed at 5 nonzero values of H for case II. The 

TABLE L Exchange parameters used for the 4 zero-field cases. 



Case 


Jo 


Jl 


J2 


I 


1.0 


-1.0 





II 


10.0 


-1.0 





III 


10.0 


-1.0 


0.1 


IV 


10.0 


-0.3 


0.001 




isotropic case I (|Jo| = \Ji\) was previously investigated 
using the MC histogram method with periodic boundary 
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FIG. 1: Results of applying the order-paramater cumulant 
crossing method to estimate the critical temperature Tn for 
case II. The straight lines correspond to linear fits to the data 
with ln'^(L7L) < 2.2. 

conditions for a,ii N — L x L x L system 0, 0, 21 1 ■ Dif- 
ferent sets of critical exponents characteristic of different 
universality classes were obtained in these previous stud- 
ies, serving to illustrate the sensitivity of analyzing MC 
simulation data for the frustrated Ising model. Cases II, 
III, and IV correspond to systems with anisotropic quasi- 
ID exchange interactions (with | Jo| ^ |Ji|) for different 
values of Jo- The anisotropic case IV was previously 
investigated|13j using the CHB algorithm and estimates 
of the critical temperature and the critical exponents /3 
and v were obtained using a data collapse method. 

The CHB algorithm employs open boundary condi- 
tions along the c-axis and periodic boundary conditions 
in the other directions. For systems with quasi- ID 
exchange interactions the CHB algorithm is more effi- 
cient than the Metropolis algorithm and allows for bet- 
ter statistics when simulating larger lattice sizes. We 
considered anisotropic lattices with N = L x L x IQL 
and L — 9, 12, 33, 36. These system sizes are smaller 
than those employed in Ref.[l3]- However, a significantly 
larger number of MC steps were used in any particular 
run. Averages were also performed over 10 independent 
simulations using different random initial spin configura- 
tions. Estimates of the errors were obtained by taking the 
standard deviation from the different simulations. Runs 
of 1 X 10^ — 8 X 10^ MC steps were used for equilibra- 
tion, and 5 X 10^ - 1.2 x 10^ MC steps were used to 
calculate averages of several thermodynamic quantities 
including the order parameter O, susceptibility xit spe- 
cific heat C, energy cumulant Ue, order-parameter cu- 
mulant Um and first logarithmic derivative of the order 
parameter Vi , as defined in Ref. [l9l | . The primary order- 
parameter O is defined in terms of the Qt^ Fourier com- 
ponent of the spin density as O = |^ ■ Ui exp(iQ • R^) | /N 
with Q — (27r/3, 27r/3, 0) in units of the lattice constants 
associated with the simple hexagonal structure. In addi- 
tion, the temperature dependence of the secondary order- 
parameter O' with the wavevector Q = (0,0,0), corre- 
sponding to the uniform magnetization, was also exam- 
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FIG. 2: Finite-size scaling of the order parameter O for case II 
where error bars on the data represent the standard deviation 
of the ten runs. 

ined. Relevance of this component of the spin density (in 
zero apphed field) on the nature of the critical behavior 
in these systems has been previously speculated ^ [l^l ■ 

Accurate estimation of the critical temperature is an 
essential first step in utilizing the histogram method for 
determining critical exponents. Temperature sweeps for 
the L = 24 system were initially performed using fewer 
MC steps to obtain a rough estimate of the transition 
temperature T/v by locating the maxima of the suscep- 
tibility xi and specific heat C. For example, for case 
IV, temperature scans from T = 6.8 — 8.0 in increments 
of AT = 0.05 yielded an estimate of Tjs = 7.40. To 
determine the critical temperature more accurately, sev- 
eral histograms were generated at temperatures above 
and below Tj^. In this case, histograms were generated 
at T =7.34, 7.36, 7.38, 7.40, 7.42 for each lattice size 
L. The order-parameter cumulant crossing method was 
used for the sets of parameters of Table I and an illus- 
trative example of the results for case II is presented in 
Fig. [T] For each case we plot the temperature at which 
Um for lattice size L' intersects with the cumulants for 
L — 9, 12 and 15. Linear fits are made using the results 
in the asymptotic region (i.e., ln~^(L'/L) < 2.2) and an 
estimate of the critical temperature T/v is obtained from 
the average value of the crossing points. The error ± ATat 
represents the standard deviation of these values. 

The fourth-order energy cumulant evaluated at T/v as 
a function of system size L was also calculated in each 
of the four cases. The results extrapolate to = 
0.666663(3) for L ^ oo in all cases, as expected for a 
continuous phase transition. Finite-size scaling analysis 
at the critical temperature T/v of the various thermody- 
namic quantities O - L-f^/", x ^ L^'"", and Vi - L^'", 
were used to obtain estimates of critical exponent ratios. 
Representative results are shown in Fig. [2] for the or- 
der parameter in case II where the error on the exponent 
ratio is obtained from the robustness of the fit. Our re- 
sults show that the secondary order parameter O' is not 
relevant at T/v for H = 0. 

In Table |TT] the zero field estimates of the critical tem- 



TABLE II: Comparison of the critical temperature and criti- 
cal exponents determined in this work with other MC studies. 







p 


7 


V 


I 

Bunker et al.[20] 
Plumer et al.[7| 
Heinonen et al.^j 


2.926(3) 
2.920(5) 
2.9298(10) 

z.oo 


0.344(7) 
0.311(4) 
0.341(4) 


1.31(3) 
1.43(3) 
1.31(3) 

1. lOl^O ) 


0.671(9) 
0.685(3) 
0.662(9) 


II 


10.651(4) 


0.355(6) 


1.33(3) 


0.670(7) 


III 


11.998(4) 


0.358(5) 


1.28(3) 


0.677(3) 


IV 

Koseki et al.[13j 


7.406(6) 
7.34(4) 


0.362(7) 
0.21(1) 


1.35(2) 
1.31(3) 


0.673(4) 
0.70(3) 


XY 




0.345 


1.316 


0.671 


Tricritical 




1/4 


1 


1/2 



perature and critical exponents determined in this work 
are compared with those obtained by other MC stud- 
ies. For each case studied, our set of critical exponents 
indicate 3D XY universality. This is in agreement with 
symmetry arguments as well as the MC results of Ref. 
but contrast with the tricritical behaviour seen in Refs. 
and where less intensive simulations were used. 
In order to estimate errors due to the uncertainty in the 
critical temperature T/v, finite-size scaling analysis was 
performed at temperatures slightly above and below T/v. 
In Fig. [21 results are shown for the critical exponents /3, 7 
and V versus the choice of critical temperature for cases I, 
II, and IV. We find that the cases with quasi-lD exchange 
interactions (cases II, and IV) are more sensitive to the 
choice of critical temperature than the case with isotropic 
exchange coupling (case I). These results indicate that 
the critical temperature region is extremely narrow for 
the quasi- ID models, highlighting the need for care in the 
analysis and interpretation of both computational and 
experimental data for such systems. 

TABLE III: Critical temperature and critical exponents for 
several points along the paramagnet phase boundary. The 
errors on the exponents are estimated from fits performed at 

r = rjv±Arjv. 



H 




/3 7 


V 





10.651(4) 


0.35(2) 1.33(5) 0.670(7) 


0.10 


10.651(5) 


0.31(3) 1.3(1) 


0.63(1) 


0.25 


10.685(5) 


0.28(3) 1.22(6) 


0.56(1) 


0.50 


10.756(5) 


0.27(4) 1.02(9) 


0.50(2) 


1.00 


10.891(5) 


0.27(7) 1.0(1) 


0.48(4) 



The effects of a nonzero applied field were examined 
here for the quasi-lD case II only. This extends the 
work of Refs. 0, where isotropic exchange was as- 
sumed and where it is argued the effect of an applied 
field is to change the symmetry to that of the 3-state 
Potts model and hence the transition should be first order 
within mean field theory. However, a more complicated 
phase diagram emerged as a result of these earlier MC 
simulations. The same number of MC steps as in the zero 
field cases were used in the present study for equilibra- 
tion and to calculate thermal averages. At field strengths 
H — 0.01, 0.25, 0.50 and 1.0 simulation results for 
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0.50 
0.45 
0.40 
on. 0.35 
0.30 
0.25 



Case p 




I 0.344(7) 

II 0.355(6) 
IV 0.362(7) 


T . 
I 


I I 




It 5^ 








1 1 ' 





-0.004 



-0.002 



0.000 
(T-TJ/T„ 



0.002 



0.004 



FIG. 3: Values of the critical exponents (3 for cases I, II, and 
IV, obtained for different choices of the critical temperature. 
Results for case III are qualitatively similar and have been 
omitted for clarity. 

the energy cumulant extrapolate to If^ — 0.666663(3), 
suggesting that the phase transitions remain continu- 
ous at these lower values. Results for the critical ex- 
ponents at five field strengths from H — Q to H — 1 
are shown in Table IIIII The error for the critical tem- 
perature is estimated from the scatter of the crossing 
points of the order-parameter cumulant data. The er- 
rors on the critical exponents are obtained from fits per- 
formed at T = Tjv ± ATn. For example, for H = 0.25, 
the fits performed at T = 10.680 yielded f3 = 0.26(1), 
7 = 1.28(4), V = 0.56(1), whereas those at T = 10.690 
yielded /3 = 0.306(6), 7 = 1.17(2), v = 0.57(1). These 
results illustrate the sensitivity of the critical exponents 
with respect to the choice of critical temperature. The 
largest source of error on the critical exponents comes 
from the uncertainty of the critical temperature. The er- 
rors from the robustness of the fits are typically smaller. 
For small values of the applied field {H = 0.1) the esti- 
mated values of the critical exponents f3 and v are slightly 
lower than the zero field cases. At intermediate field val- 
ues H = 0.5 and H = 1.0 the magnetic phase transition 
is characterized by a set of exponents that are consistent 
with tricritical mean-field values. Monte Carlo simula- 
tions results for H = 2.0 reveal a weak first order phase 
transition where = 0.666642(3). 

The results of these extensive MC simulations using 
the CHB algorithm combined with the histogram analy- 
sis technique serve to resolve long-standing questions re- 
garding both experimental and previous modeling results 
on the critical properties of the quasi- ID ISTAF where 
both XY and tricritical behavior has been reported. This 
work illustrates the difficulty in extracting reliable esti- 
mates of exponents due to the enhanced sensitivity of the 
critical region when exchange interactions are frustrated 
and anisotropic. A major focus of the present work has 
been the careful estimation of errors due to uncertainty in 
the assumed critical temperature, a feature relevant for 
the analysis of both modeled and measured data. Our 
results make clear that the classical model system yields 



XY criticality in the case of zero applied field for a vari- 
ety of assumed exchange interactions, in agreement with 
most previous MC simulations results and with most ex- 
perimental results on Ising-like ABX3 compounds such as 
S = 1 CsNiCls . Quantum effects due to SRO inherent 
in the quasi-lD 5' = 1/2 compounds are mimicked by the 
addition of a small applied field and are shown to induce 
a cross-over to tricritical, then first-order behavior. Our 
results explain recent high-resolution neutron scattering 
experiments indicating tricritical exponents in CsCoBra 
|ll| where the strength of the extra field term had pre- 
viously been estimated to be H/ Jq = 0.05 Q consistent 
with our simulations results {H = 0.5). This conclu- 
sion illustrates that quantum spin effects can control the 
experimentally accessible effective critical behaviour of 
quasi-lD ISTAF's. This result is relevant to other sys- 
tems that exhibit strong SRO before the onset of LRO. 
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